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Abstract 

We study the influence of a fast noise on Landau-Zener transitions. We demonstrate that a 
fast colored noise much weaker than the conventional white noise can produce transitions itself or 
can change substantially the Landau-Zener transition probabilities. In the limit of fast colored or 
strong white noise we derive asymptotically exact formulae for transition probabilities and study 
the time evolution of a spin coupled to the noise and a sweeping magnetic field. 

PACS numbers: 
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I. INTRODUCTION 



Landau-Zener (LZ) formula for transition probabilities at avoided crossing of two levels 
is one of a few fundamental results of non- stationary quantum mechanics. Its rather general 
character and simplicity makes it extremely suitable for versatile applications. Traditionally 
it was applied in quantum chemistry^ and in collision theory^,—. A recent treatment of 
the experiments on the quantum molecular hysteresis in nanomagnets by Wernsdorfer and 
Sessoli^*^ was a real triumph of the LZ theory. A substantial contributions to the theory of 
spin tunnelling in these molecules was made by theorists^i^. Landau-Zener formula and its 
generalizations were recently employed also for charge transport in nanostructuresiSiiiiiSiiSiii^ 
Bose-Einstein condensates^^ and quantum computingi^. 

Extensions of the LZ theory to the case of multilevel crossing are less general. Neverthe- 
less, some of them were realistic enough to justify remarkable efforts on the side of theorists 
for their analysis. Level correlations and localization in energy space were studied in^^. The 
pioneering work by Demkov and Oshero\^ treated exactly the crossing of a single level with 
a band of parallel levels. In the work^^ Hioe and Carrol solved a problem of transitions in 
a Zeeman multiplet of an arbitrary spin S in a magnetic field with a constant perpendicular 
component and a time-dependent parallel component passing through zero value. Numerous 
generalizations of these results were foun d^°i^^i^^i^^i^"^i^^i^^i^''' . A general point of view on all 
these exactly solvable models proposed by one of author^ allowed to construct an algorithm 
for series of new solvable models. Another extensions include nonlinear LZ-mode P^i^°i^^ and 
LZ problem with nonlinear sweep^^^. To apply the LZ formula and its multi-state extensions 
to real systems it is often necessary to take into account the interaction with environment. 
Such attempts were made in a series of work a^i^^i^^i^^i^^i^'''i^^i^^i^°i^'^i^^i"^^i"^^i^^ , however the 
problem was not solved completely. Kayanuma et a/ .33i34^35 j^g^yg obtained an elegant an- 
alytic result for the diagonal white noise. The non-diagonal colored noise was considered 
by Kayanumai^ for the two-level crossing without a constant coupling term. He has found 
the transition probability in the limit of infinitely short noise correlation time. His result 
was disputed by Nishino et al^. On the basis of their numerical calculations these authors 
concluded that the transitions due to non-diagonal noise are weak in many realistic situa- 
tions. The discrepancy is associated with different limiting procedures in calculations. In 
the limit of infinitely short correlations but finite amplitude of sweeping field the transition 
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probability was found to be vanishingly small. 

The influence of the noise onto the multilevel crossing was studied so far in only one work 
by Saito and Kayanumai^, who considered the 3-level crossing at a special relations between 
parameters in the limit of strong decoherence. 

The purpose of this article is to present a systematic study of the influence of noise, 
including the colored noise, onto the LZ transitions and to generalize it to multistate LZ 
problems. We demonstrate that the LZ transitions are sensitive to the colored noise much 
weaker than usual (5-like white noise. The latter can be considered as a limit of a noise 
whose correlation time goes to zero and simultaneously its square of amplitude goes to 
infinity, so that their product remains a constant. We prove that such a white non-diagonal 
noise always leads to equal population of the crossing levels. However, the noise, whose 
correlation time goes to zero, but its amplitude remains a constant, produces non-trivial 
transition probabilities as it was first found by Kayanuma-^ for a special type of the noise 
correlation function. Another subtle problem is the order of limiting processea^i. The 
resulting probabilities depend crucially on what happens first: time asymptotically goes to 
infinity or correlation time goes to zero. Analysis of these problems allowed us to reconcile 
works^^ and'^''. In our work we first find simple analytical result for a transition produced 
by a most general short-time correlated noise in 2-level systems and the change of the LZ 
probabilities produced by such a noise. We check these analytical results by numerical 
calculations. We also study the influence of the noise on transitions at multilevel crossing. 

The plan of the article is as follows. In section II we generalize the result of Kayanuma 
for transverse noise^ to the case of the arbitrary Gaussian noise in all three directions. Next 
we demonstrate its generalization to a three level system. In the section IV we study the 
time dependence of the density matrix with LZ transitions stimulated by fast noise. In the 
fifth section we propose a formula that incorporates constant transverse magnetic field and 
compare its predictions with numerical simulations. In section VI we consider the master 
equation for an arbitrary spin placed into a regular varying field and noisy magnetic field 
along the 2;-direction and a constant field along x-direction and find simple expressions for 
transition probabilities in the limit of a strong decoherence. In section VII we perform 
similar calculations for a charged particle on a periodic chain driven by a time- dependent 
electric field and compare our results with those for a completely coherent evolution. 
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II. COLORED NOISE IN TWO LEVEL LZ TRANSITIONS. 



LZ transitions in a two- level system with a non-diagonal noise were studied by Kayanuma 
in^. The Hamiltonian of the problem was chosen to be 

H = pta^ + T]^a^ (1) 

where 1]^ is the noise field with the correlation function (?7x(^i)^x(i2)) = J^e"'^"''*^"*^' and (Tj 
are Pauli matrices. In the limit of infinitely short correlation time A ^ oo Kayanuma has 
found a simple analytical result for the transition probability. 

The choice of Kayanuma corresponds to the spin 1/2 problem with noisy magnetic field 
along the x-axes only. We generalize the Kayanuma model introducing a more general 
Hamiltonian with all 3 components of random magnetic field being non-zero and with a 
most general form of the short-time correlation tensor: 

H = ptaz + '^r]i{t)ai, i = x,y,z (2) 

i 

{V^itMt2))=9^JiX\tl-t2\) (3) 

We assume that though Qij can be different for different i,j, they are of the same order 
of magnitude. We consider the limit of fast noise with A — > oo, where A is the inverse 
characteristic decay time of the correlator Qij. 

The density matrix elements for the system with the Hamiltonian obey the following 
system of ordinary differential equations: 

ip = 2(-pi2(?7^. - ir]y) + p2i{Vx + iVy)) 

ipi2 = 2{f3t + r]^)pi2 - p{r]:^ + iT]y) (4) 

ip2i = -2{(3t + r]^)p2i + p{'qx - ir]y) 
where p = pu — p22- By elimination of non-diagonal matrix elements equations (0)) are 
transformed into one integral-differential equation for the occupation difference p: 

The solution of this equation can be formally found as an infinite series in powers of 
rji that must be averaged over noise realizations. A typical term contains the product 



4 



'nh{ti)rji^{t2) ■ ■ ■ Tji^itn) ■ Its average is equal to the sum of all possible products of pair 
correlators since we assume the noise to be Gaussian. Kayanuma)2£ has shown that in the 
limit of very fast noise only the term in which the pairing is ideally ordered in time, i.e. 
the pairs are (12)(34)...(2n-l,2n), contributes a finite value into the integral. Other pairings 
contribute terms, which are by a power of infinitely small parameter 1/A smaller and can 
be neglected. This Kayanuma's observation is completely analogues to the theorem proven 
by Abrikosov and Gor'kov in their theory of impurities in a metal^^. Note that the diagonal 
component of noise is inessential in this approximation and can be omitted. These facts 
allow to write down the integral-differential equation for the average value of p as follows: 

d{p)it) 
dt 







1 COS (/ 


3{e-ti)) 





4 / cos[(3{t'-ti))F{\\t-t^\){p){t,)dt^ (6) 



where F = g^x + gyy Now we can employ the approximation of the fast noise assuming 
that the average (p) (t) almost does not change in the interval of time 1/A and that integral 
of correlation function is convergent. In this approximation we can extract {p){t) from the 
integral in the right-hand side of equation (jH)) and expand the argument of the cosine near 
the end point ti = t of the integral. The resulting differential equation reads: 

^ = -4F(2/3t)(p) (7) 

where F{q) is the cosine Fourier transform of the function F: 

F{q) = / cos(gr)F(A|r|)rfr (8) 

J —oo 

Note that the characteristic value of q are of the order A and, respectively, essential values 
of t in equation ((7j) are t ~ A//5. We see that essential values of time go to infinity together 
with A. It shows that the order of limiting processes is indeed very important. Here we 
first calculate the transition probability, i.e. the diagonal elements of the density matrix at 
t = oo and at very large, but still finite A and put it infinity in the end. Solving equation 
dZj), we find: 

(p)(t) = (p)(-oo)exp (-4 f_^F{2(3t')dt'^ (9) 

Together with the standard equation pn + P22 = 1 equation © determines average occupa- 
tion number of each level at any time. The most interesting are the transition probabilities 
at t = +00, which can be obtained from the same equation 0. Note that the integral in 
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the exponent at t = +00 becomes equal to 

F{2(3t')dt' 



TT 

2P 



m 



(10) 



Thus, from equations Q and (fTUj) we find: 



Pii(+oo) 



/)22( + 00) = - 



1 + (2pn(-oo) - l)exp 
1 - (2pii(-cx)) - l)exp 



Putting pii(— 00) = 1, we find the transition probabihty: 



Pi- 



P22( + 00) |pii(-oo)=l= 77 



1 — exp 



27rF(0) \ 
2nF{0) \ 



27rF(0)Y 



(11) 
(12) 



(13) 



For the standard white noise the correlators {r]i{t)r]j{t')) turn into delta-functions in the 
limit A = 00. It means that their values at t = t' become infinitely large. In particular 
it means that F{0) is infinitely large for the standard white noise. In this it is 

seen from equations (fTT|) . (fT^ the occupancies of both levels are equal to 1/2. Thus, the 
standard white noise leads to complete loss of initial state memory and equipartition of the 
levels. On the other hand, if the amplitude of noise remains finite, the occupation numbers 
conserve memory on the initial state. 

In the limit of the fast noise the transition probabilities are determined only by the average 
square of non-diagonal noise and are not sensitive to the diagonal noise. To illustrate this 
statement we consider a special case when the diagonal noise does not correlate with the 
non-diagonal one. Then, as it is seen from equation (0), the averaging over the z-component 
of the noise leads to multiplication of the coefficients in this integral-differential equation by 
the Debye- Waller factor 



(exp ^ r/^(r)c/rj) 

= exp (^-j^{t - ti) J2^gU0)d9 

An essential interval of integration over ti near ti = t is 1/A. In this interval the Debye- 
Waller factor changes by the value ~ 1/(A)^ and with this precision is equal to 1. 

For a special case gxx{T) = 9yyi''~) = we reproduce the Kayanuma's 

result*^. 
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III. NONDIAGONAL NOISE IN SPIN-1 LZ THEORY. 



The Hamiltonian if of a general multi-state LZ problem (see for example^S ) has the 
following matrix form 

H = Bt + A, (14) 

where 5 is a diagonal matrix and matrices A and B are independent on time. However, 
in a situation of a general position only two levels cross. Several levels can intersect at the 
same moment of time only due to a special symmetry. Such a symmetry is systematically 
realized in the model of an arbitrary spin S ^ 1/2 placed into an external magnetic field 
that has a time-dependent z-component vanishing at some moment of time t = and a 
constant transverse componentAS. The corresponding Hamiltonian is: 

H = ptS, + gS^ (15) 

where jS and g are constants. This exactly solvable model for a spin higher than 1/2 was 
employed in the theory of Stark effecli^; some recent applications can be found in^S*^. 

In this section we generalize the result to a spin 1 system. We consider a following 
Hamiltonian for the spin 1 system in a random magnetic field: 

H = {(3t)S, + r],S, (16) 

i=-l 

where Si are the spin projection operators for S = 1. The density matrix p depends on 8 
independent parameters. We denote p+ = pi^i — po,o and p_ = P-i,-i — Po,o- The derivation 
of equations for p+ and p_ can be done in the same spirit as it was done for the spin 1/2 
case. Looking for a solution of the evolution equation for density matrix in the form of 
perturbation series and retaining after the averaging over the noise only the leading terms 
in 1/A, we arrive at an integral-differential equation. Additional care must be paid to the 
integrals of the form j'i^dti j'^^dt2exp{i(3{tl + tl))gij{X\ti — t2\) that did not appear in 
the two level system, but appeared in the series for the spin S = 1. One can check that 
this integral is of the order 1/A and hence we disregard it and all terms that contain it. 
After lengthy but straightforward calculations we find that in the leading order in 1/A the 
elements p+(t) and P-(t) satisfy the following integral-differential equation: 

-2t^dt,cos{(3{t^-tj)/2)F{X\t-U\) ' \ \ ^"^"^^ I (17) 





The approximation of the fast noise allows to transform this equation into the following 
differential one: 



d 
~dt 



-it) i 



-2F{(3t) 




P-ih 



(18) 



Thus, the problem is reduced to a linear differential equation with a constant matrix coeffi- 
cient. The formal answer is: 




exp 



if F{(3\T\)d7 
J —oo 













^P-(-oo) J 



(19) 



This formal solution can be transformed to a more explicit form with the help of the matrix 
identity: 

^11/2 



exp 



-7 



1/2 1 



^ cosh7/2 — sinh7/2'^ 



V 



sinh7/2 —cosh 7/2 



(20) 



The transition probabilities are defined by this matrix acting on vectors 



if initially 



the projection +1 was occupied , or 
probabilities are: 



if the projection —1 was occupied. The transition 



Po^o = i (1 + 

3 



1 + ie-37/2 I ie-7/2 
3 -r gts -r 



1 + ip-37/2 _ ip-7/2 
3 ^ 6^ 2^ 



(21) 



Po^+ - Po^- - I (1 - 
where 7 = 7rF(0)//3. 

We see that in the case of 3 levels the main results obtained in the previous section for 
a 2-level system persist: the standard white (7 = 00) noise leads to equal population of 
all 3 levels, whereas the fast noise with a finite amplitude results in non-trivial transition 
probabilities. 



IV. TRANSITION TIME FOR COLORED NOISE 



In section 2 we have shown that the typical time for establishing the asymptotic of the 
transition probabilities is A//3F(0). It is useful to look at the same problem from a different 
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point of view. Namely, we will analyze the behavior of transitions driven by the standard 
5-like white noise, their typical rates and times. If the action of the standard white noise is 
limited by some finite time interval, it becomes physically equivalent to the fast noise with 
a finite amplitude. Simultaneously we will study directly the influence of the white noise 
onto the LZ transitions. 

Let us consider the problem with the very beginning. The Hamiltonian is the generator 
of a random rotation: 

H = h(t)S (22) 
where h{t) = ho(t) + i]{t) and ho(t) = (3tz + Tx. Here we consider the isotropic white noise 

The density matrix as any Hermitian S x S matrix can be represented as a sum: 



+ ... (23) 



where the last term in equation (j29|) contains 5* operator factors. All irreducible tensors 
p^^^; j = 0...S evolve separately. The scalar p(°) = {l/N)Trp is a constant. The vector p(^) 
obeys an obvious equation: 

p(i) = -h X p(i) (24) 

The modulus of the vector p<i) is conserved. 

Dynamic equation for the 2-nd order symmetric tensor pik reads: 

Pik ^imnhmPnk ^kmnhmPin (^^) 



The extension for the rest of irreducible components is obvious. In equation (j24j) for the 
vectorial part it is useful to apply the interaction representation: p^^\t) = Uo(t, to)p*-"'"^(t), 
where Uo{t,to) is the evolution matrix for the field ho(t) = z^t + xT. The reduced vector 
p^^^ (t) obeys the following equation: 

|p«(t) = -r^xp~«(t) (26) 

Here fj = Uq^^t] x Uq). In terms of the Cartesian coordinates this equation reads: 

Pi = TikipkVi (27) 

where the matrix Tiki is: 

Tikl{t,to) = {Uo)mi{Uo) nk^mln 

(28) 



Initial values of p and p coincide since Uo(to, to) = I- Equation (jTfjl can be solved by power 
expansion over the noise: 

P^it) = piito) + f Tiki{t' ,to)rii{t')dtpk{to) + ... (29) 

When averaging the expansion over the noise, all odd terms vanish. To understand what 
happens with even terms consider first the quadratic term: 

I dti I dt2Tik^i^{ti,to)Tk^^k,i2{t2,to){Vh{ti)vi2it2))pk{to) 

Jto Jto 

For the isotropic noise the contribution to the vector p{t) up to the second order in rj can 
be represented as follows (we write equations for averages omitting the angular brackets): 

rt 
'to 



P2,i(t) = {Sik + dhMikih)^ Pkito) (30) 



where 

Mikit) = {^/2)Timlit, to)Tmklit, to) = -^5ik (31) 

Here we have used properties of the orthogonal matrices Uq: (f/o)mj(t^o)mfc = ^ik- 
Next we proceed to the quartic term: 

rt rti rt2 rtz 

pA,i = / dti / dt2 / dt^ / dt4Tipk(ti,to)Tpgi(t2,to)Tgrmi't3:'to)Trsnit4:,to) 

J to J tQ J tQ J £() 

X {'nk{tl)Vlit2)Vrn(t3)Vn{U))Ps{to) (32) 

Quartic average of Gaussian field decays into quadratic averages: 

{'nk{ti)rii{t2)nm{tz)rin{ti)) = {r]kiti)r]i{t2)) {r]rrXt3)VniU)) 

+ {Vk{tl)r]m{t3)){r]i{t2)r]n{t4^)) + {r]k(ti)r]n{t4)){r]i{t2)r]m{t3)) (33) 



Only the first term contributes to the integral ()32|) . Two others are zero at correct time 
ordering ti > t2 > t^ > t^. In both cases considered earlier we obtain: 



h,i= I dti f dt2Mij{ti)M,k{t2)Pk{to) (34) 

Jto Jto 

The same result could be obtained from an effective equation of motion for p: 

Pi = MikPk (35) 
Thus, the operator M plays the role of effective non-Hermitian Hamiltonian. 
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For the white isotropic noise the solution of the equation ()35|) is: 



Pi(t) = Pi(to)exp{--ft) 



(36) 



At t — s> oo the formula ()36p always leads to the occupation numbers p = 1/2. As we know, 
this does not happen for the colored noise with a finite amplitude. The reason is that in the 
genuine LZ problem the solution strongly oscillates with the frequency roughly u!{t) ~ f3t 
long before and after the level crossing point. This introduces a new energy scale that must 
be compared with A. For time in the range \(3t\ < A, the approximation of white noise is 
roughly valid even for finite amplitude noise, but beyond this interval of time the oscillations 
of the LZ solution become faster than the correlation time of the noise, and the action of 
the noise is suppressed by the oscillations. 

In Fig.l we demonstrate a typical evolution of the probability for the system to stay in the 
same state as a function of time. The evolution reminds diffusive motion that slowly stops 
at large absolute values of time. The sharp change of the absolute value of the amplitude 
near t = is due to the constant transverse field. 

To estimate roughly the transition probability one can apply the standard white noise 
approximation in the time interval \t\ < CX/jS and accept that at \t\ > C\/(3 no transitions 
due to the noise happen. The parameter C is a constant of the order of unity. Then we 
automatically get the result p{oo) = exp(— 2CA7//3)p(— cxd). According to the definition, 
7 ^ F{0)/X in agreement with the calculations of the previous sections. Summarizing, the 
transitions mediated by the fast noise proceed during the large time interval of the order 
of A//3. This is the reason why its total effect remains finite, though it is very small on a 
typical time scale of the usual LZ-transitions mediated by constant field (r ~ r//3)^. 

V. LANDAU-ZENER TRANSITIONS IN A CONSTANT TRANSVERSE FIELD 
AND IN THE COLORED NOISE. 

In previous sections we assumed that only the noise was responsible for transitions, 
whereas the regular part of the Hamiltonian operator was diagonal. In this section we 
incorporate a regular non-diagonal operator (the transverse field) into the Hamiltonian to- 
gether with a noise. The most general Hamiltonian for such a two-level system reads: 



H = pta^ + gux + '^r]i(t)ai, i = x,y,z 



(37) 
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FIG. 1: Typical evolution of the probability P[t) = |ai(t)p, where ai{t) is the amplitude of the 
first diabatic state. The choice of parameters is as follows, gij = 5ijJ'^e~^^^~^'^, /3 = 1, A = 80, 
J2 = 0.18, r = 0.7. 

where rji are coordinates of the vectorial Gaussian noise with correlators given by (jHl)- 

The second term in the Hamiltonian ()H7|1 can be considered as a constant transverse 
magnetic field acting on a spin S = 1/2. The Hamiltonian (j37|) describes a spin (q-bit) 
weakly interacting with the environment, for example with the nuclear spin bath3S. If this 
interaction is so small that the bath relaxation is much faster than the inverse interaction 
energy, the bath can be treated as a fast noisy magnetic field. The measurements of the 
LZ transition probabilities can provide an information about the strength of the coupling to 
the bath. Another possible example is a molecular nanomagnet in fluctuating dipole field. 
The Hamiltonian (jHTf) may be relevant to the quantum shuttle problem where avoided level 
crossings occurred to be important^i. The fast noise in this example corresponds to thermal 
fluctuations. 
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As it was discussed in the previous section, in the interval of time \pt\ « A the spin 
experiences an equivalent white noise with the amplitude 7 = F{0) / (2A) vanishing at A — *^ 
00. The noise causes a slow decay of the population difference p = pu — p22- The action 
of the noise becomes noticeable only after long evolution time of the order |A//3|. On the 
contrary, the transitions due to the constant transverse field proceed predominantly during 
a much shorter interval of time \t\ < \g/P\. During this interval of time the role of noise is 
negligible comparing with the role of the constant transverse field. Let us choose to much 
less than | A//3 1 , but much larger than g/(3 to make sure that the transverse magnetic field is 
neglegible beyond the interval of time (—to, h)- The evolution from —00 to —to is influenced 
mainly by the noise. The result of such an evolution is almost the same as that in the absence 
of the constant transverse field with to = 0. Since the effective attenuation coefficient for p 
is an even function of time (see equation ((7j)), the resulting exponent for p(to) is twice less 
than that for the problem of the section (2) at t = +00 (equations ()9I10|) : 

p{-to) = e--^(°)/^ (38) 

In the interval of time |t| < to the noise is neglegible and transitions are completely 
due to the genuine Landau-Zener mechanism, i.e. due to the constant transverse field g. 
Since to ^ g/P, we can use the LZ formula to determine p(to)- Let ai(t) and a2(t) be 
the amplitudes to find the system at the first and the second diabatic states, respectively, 
(to find the spin projection equals 1/2 or —1/2) at the moment of time t. Their values at 
t = —to can be expressed in terms of the population difference p* and two phase factors as 
follows: 

a,{-to) = e^'^^^{l+p*)/2 ^^^^ 

a2(-to) = e^<^y(l-p*)/2 
The phase diference (pi — 02 is essentially random and independent on p* since, the solution 
of the Shrodinger equation strongly oscillates at t < —to and transitions are presumably due 
to the non-diagonal noise. The amplitudes at time to are related to those at t = —to by a 
linear relation aj(to) = Sij{to, —to)aj{—to), where 5* is the transition matrix in the noiseless 
LZ Hamiltonian from t = —to to t = to- Using this property and averaging over the random 
phases 0i and 02 and over p* with < p* >= p{—to), we arrive at a following result: 

p{to) =< \ai{to)\' - \a2{to)\' >= e-'^^W/'3(2e--l^l'/^ - 1) (40) 
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FIG. 2: The probability to stay on the same diabatic state in a two level system at constant 
coupling to the noise as function of transverse magnetic field. The choice of parameters is gij = 
5*x5jx.j2e-^l*-*'l, = 0.28, A = 125, P = 1. 

Here we employed the fact that our choise of time to allows the approximations 
\SijitQ, —to) I ~ \Sij{oo, —oo) \ and p(— to) ~ p(0) The evolution from to to t = +oo brings an 
additional exponential factor equal to that in equation (j38|l : 

p(t ^ oo) = e-2-^W//^(2e-'^l^l'/^ - 1) (41) 
Correspondingly the probability to stay on the same diabatic level is 

Pi^i = i(l + e-2-^(o)//3(2e--lfl'/'3 - 1)) (42) 

In the limit of the noise only [g = 0) and of zero noise {F[0) = 0) the result ()42|) transits 
into formula ()13|) or to the LZ formula, respectively. 

The matching procedure used in this section is asymptotically exact at A ^ oo. To 
check how it works at large but finite A we studied the LZ transitions subject to a fast noise 
numerically. 

We simulated the evolution generated by the Hamiltonian (pTTj) . For simplification we put 
Vy = Vz = and accept the following form for the correlator of the noise x-component: 
{Vx{t)rix{t') = Jjexp(— A|t — t'|). The time interval of the evolution was chosen to be much 
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larger than \/(3. In Fig. 2 we depict the probabihty to stay in the same state after the 
level crossing vs. the constant transverse field at a fixed coupling to the noise Jx = 0.28. 
Each discrete point represents the averaging over 100 simulations with the same coupling 
constants. The solid line is the graph of the theoretical formula (j42|) . The deviations of 
the simulation results from the analytical predictions do not exceed the accuracy of our 
calculations. We conclude that equation ()42|) describes well the transition probability for 
the LZ system subject to a fast noise. 



VI. ARBITRARY SPIN IN A STRONG DIAGONAL RANDOM FIELD 

We have demonstrated that the non-diagonal 5-like white noise leads to the equilibration 
of population between all states. It is not correct for the white noise directed along the 
sweeping field. Such a noise does not couple different states and its action leads to the loss 
of coherence only. As it was shown earlier, in the case of a 2-level system it results in a 
Debye- Waller factor for p. We consider the case of a general spin S placed into a regular 
field ho = zl3t + xT and the random field directed along 2;-axis. Its Hamiltonian reads: 

Ht,t = H + ff„„,,, = ptS, + TSx + r/(t)5, (43) 

We assume < ■r]{t)r]{t') >= 2^5{t - t'). 

Following2^, we expand the solution of equation ()43p in the power series over the noise 
amplitude and average each term. The resulting series over powers of 7 is a formal solution 
of a differential equation known as master equation?^ 

^ = -i\H{t),p{t)]-^[S^,[S-,,p{t)]] = 
-i[vtS, + r^,, p{t)] - i{Slp{t) + p{t)Sl - 2S,pS,) (44) 

It is convenient to introduce notations Tij = T < i\Sx\j >■ Below we write down equation 
(j44p for diagonal and non-diagonal elements of the density matrix separately: 

Pa = —i{^i,i-l{Pi-l,i — Pi,i-\) + ^ i,i+l{Pi+l,i ~ Pi,i+l)) 

i = —S, —S + 1, . . . , S 



(45) 



Pmn = -{iPt{m -n)+ 7(n - m)2)p„„- 

^{^^ m,m—lPm—l,n Pm,n~\-1^ n+l,n ~l~ ^ m,m-\-lPm+l,n Pm,n—l^n—l,n)j ^ 7^ ^ 



(46) 
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It is possible to find an asymptotically exact solution of equations (j^ . PU]) in the limit of 
strong noise 7 ^ F, /9. In this limiting case, non-diagonal elements of the density matrix are 
~ r/7 times smaller than diagonal ones. Indeed, let us disregard the dynamical term pij 
in equations for the non-diagonal elements. We will justify this approximation later. Then 
equation (|46|) at n = m ± 1 implies that the non-diagonal elements Pn,n±i are suppressed 
comparing with diagonal matrix element by the factor ~ r/7. The matrix elements Pn,n±2 
are suppressed by the same factor with respect to pn,n±i etc. The characteristic time interval 
following from equation (jl^jl is At ~ \pnn/ (XPn,n±i) ~ From this estimate we find that 

the time derivative of the largest non-diagonal matrix element |pi,i±i| ~ (r^/7)|pj,j±i| 
1Pi,i±i can be neglected. Retaining only main diagonal and two adjacent non-diagonals in 
the matrix equations (j45j) . (j46|) . we express the non-diagonal elements in terms of diagonal: 

Pi+l,i = ~ i/3f+^' iPi,i ~ Pi+l,i+l)) 

Pi,i+l = ''fi^t^^ {Pi+l,i+l ~ Pi,i)^ i^'^) 

Pij = 0, > 1 

The problem is reduced to determining of the 25* + 1 diagonal elements. Let introduce a 
vector c with the coordinates Cj = pa. They are probabilities to find the spin in a particular 
eigenstate of the operator S^- Substitution of (jlTj) into (jlSj) gives a differential equation for 
the vector c{t) 

c{t) = (— ^ — ^ )Mc{t) (48) 

were the constant matrix M has following matrix elements: 



Ma = -{<i + l\S^\i >2 + < i\S^\i - 1 >2) = -1(^2 ^s-e) 



(49) 



Mi,i+i = Mi+i,, =< i + l\S^\i >2= \{S + i + 1){S - i) 
All other elements are zeros. Equation ()48|) can be easily integrated 

To find transition probabilities we take limits of integral in equation as to = ~ 00 and 
t = +00: 

27rr2 

c(+oo) = exp(— — M)c(-oo) (51) 

P 

This result demonstrates that, for a general spin S as well as for a spin 1/2, the transition 
probabilities do not depend on the specific value of 7 provided that 7 is large. Below we 
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present explicitly the transition probabilities for some values of spin. We denote Ei = 
g-r^TT//?^ ^2 = e-^^''^/^ and E3 = e-^^^^lf^. Then for 5 = 1/2 the formula (HH) reads 

A/2^1/2 = ^-1/2^-1/2 = 1(1 + ^1) ^^2) 
-Pi/2^-1/2 = -P-1/2^1/2 = |(1 — E\) 

For 5 = 1 

-Pi^o = -P-i^o = -Po^i = -Po-»-i = |(1 ^ E'l) 
Po_o = |(l + 2^2) 

The results fl52|) and ()53|) coincide with already known solutions for two- and three-level LZ 
models with strong decoherence^. A new result for a spin S = 3/2 that follows from (|51|1 
reads: 



P3/2- 


y3/2 - 


-- P-3/2^ 


-3/2 — 


i + ^^3 + i^2 + 


20^1 










P3/2- 


.1/2 = 


-- P-3/2^ 


-1/2 = 


-Pl/2^3/2 = P-1/2^ 


-3/2 = 


1 
4 




\E2 + 
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Pl/2- 


.-1/2 


= P-1/2- 


-.1/2 = 


\ ~ 20^3 + \E2 - 













(54) 



In Fig3 and Fig A we show the dependences of transition probabilities for 5 = 3/2 on F. 
In adiabatic limit / (3 » 1 all states are equally populated after the evolution. 

VII. ELECTRON MOTION DRIVEN BY ELECTRIC FIELD 



Now we consider another multi-state Landau-Zener model that, in the absence of the 
noise, was solved exactly2£. Physically the model describes the transport of a charged 
particle in a regular linear chain driven by a time-dependent homogeneous external field. 
Such a model is an idealization of atomic scale molecular wires or linear arrays of quantum 
dots. An important assumption in the model that makes it exactly solvable is that all sites 
of the chain are identical and equidistant. External electric field splits the energy levels at 
different sites of the chain and suppresses the tunnelling between them. Hence the transitions 
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'"'3/2- 






->1/2 


P 




3/2- 


->-1/2 


p 




3/2- 


->-3/2 




FIG. 3: . 

Transition probabilities from the state with = 3/2 to any other state as functions of transverse 

magnetic field. 

proceed in a narrow time intervals close to moments at which the electric field becomes zero. 
The noise in such a system arises due to thermal fiuctuations chaotically changing the energy 
of the electron. We suppose that there are no correlations of noise at different sites. 

Let us denote \n) a state located at the n-th site of the chain. We assume that these 
states form a complete orthonormal set (Wannier basis). In terms of this set the electron 
Hamiltonian with linear dependence of an external field on time reads: 



H — I + 1 I +C.C.) + nvt \n){n\ +r]n{t) \ n){n 



(55) 
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FIG. 4: Transition probabilities from the state with Sz = 1/2 to any other state as functions of 
transverse magnetic field. 

where v and g are constants and we assume that the noise power is the same for all sites 
i.e. < Timlin >— '^l^mn- The derivation of the master equation for this case is similar to 
that of the previous section. We consider only the limit of strong noise j » g. Then, as in 
previous example, the non-diagonal elements of the density matrix pi j with |i — j | > 1 can 
be neglected. Equations for diagonal matrix elements of the density matrix are: 

Pnn — —ig{Pn+l,n + Pn-l,n — Pn,n+l — Pn,n-l) (56) 



1/2->3/2 



1/2->1/2 



1/2->-1/2 



1/2->-3/2 



19 



Equations for non-diagonal elements after the averaging over the random noise read: 

Pn+l,n = {-ivt - 27)p„+i - ig{pn,n - Pn+l,n+l) 

(57j 

Pn-l,n = {iVt - 27)p„,_i,„ - «fi'(p„,„ - p„_l,„_l) 

Neglecting again time derivatives in these equations, we find: 

(58) 



Pn+l,n — Pn+l,n+l) 



Pn--l,n TOt— 27 (^"i" Pn— l.n— 1) 

Substituting (jKH|) into the equations (jKHjl for diagonal elements, we obtain the evolution 
equations for diagonal matrix elements of the density matrix pn,n'- 

^"•'^ = [wt^ ~ wtT2^) ~ ~ ^"-^'"^ ^^^^ 

Without loss of generality we can assume that initially, at t = —00, the particle was located 
at the site number zero. It can be treated as initial conditions for the master equation: 
Po,o{t = —00) = 1 and all other elements of the density matrix are zeros at t = —00. Then 
diagonal matrix elements pn,n acquire the meaning of transition probabilities from zeroth to 
the n-th site at a current time t. As in the previous example we can find the solution for a 
chain of arbitrary number of sites in the matrix form. 

In the limit of infinite number of sites a compact solution can be found by employing the 
Fourier-transformation p„ „ = J^"^ e^^'^u{(j),t). The system of coupled differential equa- 
tions ()59p is diagonalized by this transformation. Corresponding differential equation of the 
first order for the function u{(j),t) is readily solved. Its solution with the initial condition 
uid), —00 = 1) is 



— / — h arctan — ( 1 — cos o) 
I 2 27]^ 



(60) 



u{(f), t) = exp 

In the limit t — ^ +00 it approaches its limiting value: 

u(0,t^+oo) = e^(^-™^'^) (61) 
By the inverse Fourier-transformation we find the diagonal elements of the density matrix: 

Pn,nit +00) =Pn = e-^4(^) (62) 

V 

Here P„ is the transition probability from the site with the index to the cite with the index 
n . 
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FIG. 5: Comparison of transition probabilities in coherent (without noise) and incoherent (with 
strong noise) LZ models on the chain. 

It is interesting to compare results of this calculation which incorporates a strong noise 
with the transition probabilities without noise. In the absence of the noise (7 = 0) the 
transition probabilities are2&: 



Fig.h shows a typical behavior of the transition probabilities for both cases. The dif- 
ference in the behavior is clearly pronounced. In the absence of the noise the transition 
probabilities oscillate as functions of n and g (see^^ for details). These oscillations arise 
due to the interference among the amplitudes of different Feynman paths leading from the 
initial to the final point. In the case of a strong noise these oscillation are suppressed by 
the noise imposed decoherence and the probability distribution is a smooth bell-like curve. 




(63) 
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A simple parameter that is related to the effective diffusion coefficient and can be measured 
experimentally is the average square displacement of the particle during one sweep of the 
external field. For a chain with a strong noise it is: 



Despite a strong difference in the distribution functions, the average square displacement 
(j64|l coincides with that for the coherent evolution without noise. 

VIII. CONCLUSION 

In conclusion, we have derived the formula for the transition probabilities at the non- 
adiabatic crossing of two levels coupled by a fast non-diagonal random field in the Landau- 
Zener approximation. Depending on the strength of parameters it interpolates between 
the Landau- Zener formula for noiseless system and the Kayanuma's result for transitions 
mediated by the transverse noise only. We have determined the time intervals during which 
transitions are substantial and showed that time of transitions mediated by the constant 
field is much shorter than the time necessary for transition caused by the fast noise. Our 
numerical simulations are in a good agreement with the analytical formulas . We have 
discovered an important property of the non-adiabatic tunnelling process: its probability 
depends not on the "power" of the noise but rather on the coupling to the noise 

g{0). The former usually is responsible for the decoherence rate in systems without time- 
dependent fields. Measurements of the LZ transition probability can provide the information 
about the value of this coupling. The multi-level systems placed in regular time-dependent 
fields feel subtle differences of the noise statistical properties. 

We assumed that the diagonal noise does not dominate. In this situation it does not 
play any role. However, if the component of the diagonal noise is much stronger than non- 
diagonal ones, the situation may change drastically. The problem of strong transverse field 
compatible with noise inverse correlation time also remains open. 

We have extended the formulas and methods employed for two-level systems to a couple 
of multi-state LZ systems: an arbitrary spin experiencing the time-dependent regular and 
random magnetic field and a linear chain of sites in the external time-dependent homo- 
geneous electric field plus noise. The Landau-Zener transitions for spins higher than 1/2 




V 



V 



(64) 



n=— oo 
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were observed in a number of system a i i . Therefore we believe that our solutions can be 
checked experimentally. 
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